Non-Equilibrium Dynamics and Superfluid Ring Excitations 
in Binary Bose-Einstein Condensates 
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We revisit a classic study [D. S. Hall et al., Phys. Rev. Lett. 81, 1539 (1998)] of interpenetrating 
Bose-Einstein condensates in the hyperfine states \F = l,m/ = —1) = |1) and \F — 2,m/ — +1) = 
|2) of *^Rb and observe striking new non-equilibrium component separation dynamics in the form 
of oscillating ring-like structures. The process of component separation is not significantly damped, 
a finding that also contrasts sharply with earlier experimental work, allowing a clean first look at 
a collective excitation of a binary superfiuid. We further demonstrate extraordinary quantitative 
agreement between theoretical and experimental results using a multi-component mean-field model 
with key additional features: the inclusion of atomic losses and the careful characterization of trap 
potentials (at the level of a fraction of a percent). 

PACS numbers; 03.75.-b, 03.75.Dg, 03.75.Mn 



Binary mixtures of Bose-Einstein condensates 
(BECs) [H, i, [3 display rich phase separation behavior 
that is driven primarily by the nonlinear interactions 
between the different atomic species or states that make 
up the condensates. This places them at a remarkable 
interface between atomic physics, nonlinear and wave 
physics, and non-equilibrium statistical physics. The 
first experiment to directly probe the dynamics of binary 
BECs found complex motion that tended to preserve 
the total density but quickly damped to a stationary 
state with non-negligible component overlap. Other 
experiments have since explored various regimes. For 
example, partially condensed gases of ^^Rb exhibited 
dramatic transient non-diffusive spin polarizations 
reminiscent of longitudinal spin waves Rotating 
binary BECs were observed to evolved from an initially 
coincident triangular vortex lattice through a turbulent 
regime into an interlaced square vortex lattice In 
optical traps, BECs consisting of atoms in different 
Zeeman levels of ^^Na formed striated magnetic do- 
mains 0] between which the spin could tunnel 0] . In 
^^Rb, different Zeeman levels in a quenched condensate 
formed ferromagnetic domains of variable size with 
interesting topological features 0]. 

In understandin g c omponent separation, theoretical 
investigations [l^, [ill . [l3 . [isj have demonstrated that 
the stationary state of a BEC mixture depends criti- 
cally on the intra- and inter-species scattering lengths, 
as does its stability against excitations. On the other 
hand, theoretical studies of static and dynamic proper- 
ties [l3, 0j EBj [Si j such as the excitation spectrum 
and the nature of low-frequency simultaneous collective 
excitations, remain unverified experimentally. In this pa- 
per we begin to fill in the gaps between theory, numerical 



simulation, and experiment. 

Our experimental approach follows closely that of the 
first interpenetrating condensate experiment at JILA 0] . 
We begin with a single condensate of ^^Rb in the |1) 
state in a time-averaged orbiting potential (TOP) mag- 
netic trap [Toj . The magnetic quadrupole gradient of 
144(1) G/cm (axial) and rotating magnetic bias field of 
8.32(2) G were adjusted to confine both hyperfine states 
with coincident potential minima Q. We used RF spec- 
troscopy to establish that the magnitude of the rotat- 
ing bias field varies by no more than ±0.3% over a 
full rotation — thereby controlling one of the dominant 
contributions to the cylindrical asymmetry of the trap. 
The radial [axial] trapping frequencies were measured 
to be 30.832(13) Hz [85.90(2) Hz] for the |1) state and 
31.076(13) Hz [86.22(2) Hz] for the |2) state — a differ- 
ence of about 0.8% [0.4%] that agrees with a calculation 
that includes the effect of quadratic Zeeman terms. 

A 0.930 ms two-photon (z/^ = 4.62020 MHz, z/^w = 
6.83006800 GHz) 7r/2-pulse ^ with approximately 
500 kHz intermediate-state detuning then drives 50(2)% 
of the atoms from the |1) state to the |2) state. After the 
pulse ends, the system is allowed to evolve for a time t. 
During this evolution, a weak radiofrequency "shield" ap- 
plied at 6.3 MHz reduces the influence of trapped thermal 
atoms produced by inelastic losses. The trap is turned off 
by reducing the magnetic quadrupole to zero in ^ 200 /j,s; 
the BEC expands as it falls. Finally, 21.5 ms after releas- 
ing the atoms from the trap, we reduce the bias field to 
1.0 G and image 2l| the atoms in either state with one 
of two cameras, the first oriented to view the atoms from 
the top (along the direction of gravity and axis of sym- 
metry), and the second along a perpendicular axis. The 
entire procedure is repeated for various evolution times 
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FIG. 1: Top view of a time-sequence of experimental and numerical density profiles for — 3.50(5) x 10^ *^Rb atoms with 
equal populations in the |1) and |2) states. The first row shows the measured density profiles for the |1) atoms, while the second 
and third rows give numerical results including losses and different trap frequencies (Num. A) and without those additional 
model features (Num. B). A similar arrangement is given for the |2) atoms in the fourth, fifth, and sixth rows. The field of 
view in all pictures is approximately 100 ^m on a side. The evolution time (in ms) for each column is indicated in the top row. 



<, each with the same number of atoms N, resulting in a 
series of post-expansion top- and side-view images depict- 
ing the cvohition of each component of the condensate. 

After the transfer, the system can be described the- 
oretically by a pair of coupled nonlinear partial dif- 
ferential equations of the Gross-Pitaevskii (GP) type 
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j/m are the coupling constants for 



elastic interactions between atoms in states \i) and \j) 
with scattering lengths Cij, and Vi represents the trap- 
ping potential for species The key features that dis- 
tinguish the present model from earlier ones are (i) the 
different trap frequencies for different species, and (ii) 
the explicit inclusion and accurate quantitative charac- 
terization of the loss rates Ti. Here, Fi = -1(7111 IV"! I ''^ + 
li2\i>2?) and F2 = |(7i2|^iP+722|V'2p) represent losses 
for species |1) and |2), respectively. Physically, three- 
body recombination losses dominate in a pure 11) BEG 
(decay constant 7111 = 5.4(11) x 10"^" cm^/s [H]), but 
two-body inelastic collisions dominate in mixed and pure 
|2) BECs (decay constants 712 = 0.780(19) x 10"" cmVs 
and 722 = 1.194(19) x 10"^^ cm'^/s, values we report here 
for the first time, as obtained from loss data with mixed 
and pure condensates, respectively) . Other loss processes 
are of limited experimental significance in this particular 
system and thus are not included in the GP equations. 

The binary condensate begins in a highly non- 
equilibrium state as a result of the transfer of atoms and 
small differences in the inter- and intra-species scatter- 
ing lengths. One would then naturally expect the sys- 



tem to evolve in the direction of its stationary ground 
state, which consists of a ball of |2) atoms surrounded 
by a shell of |1) atoms [i^. Such evolution was observed 
in Rcf. 01, in which the system apparently arrived at its 
ground state without exhibiting any oscillatory behavior. 
In our experiment the motion is not rapidly damped 2J] , 
in marked contrast to the earlier results [2^. A typi- 
cal time sequence is shown in Fig. [1] which we present 
alongside the results of the corresponding numerical sim- 
ulations [2^ with and without the key model additions. 
The superb agreement of the improved model (Num. A) 
with the present experiment emphasizes one of our main 
results: the motional damping observed in Ref. [2] is not 
intrinsic to the component separation process. 

The shortcomings of the unimproved model (Num. B) 
are clearly visible in Fig. [U as it exhibits ring struc- 
tures when the experiment does not (e.g., at 60 ms) and 
does not exhibit rings when the experiment does (e.g., at 
120 ms). We find the modified timing of these features in 
the improved model is due primarily to the different trap 
frequencies experienced by states |1) and |2). Both of the 
key model enhancements also contribute to more subtle 
changes in the coetaneous atomic density distributions. 

The ring patterns in state |1) at t = 30 and 120 ms, 
as well as that of state |2) at t = 70ms, constitute the 
most striking manifestation of the dynamical evolution 
and phase separation of the hyperfine state components 
of our system. Taking one period to be the time be- 
tween the two ring-like patterns appearing at 30 ms and 
120 ms, we determine an oscillation frequency v 1\ Hz 
— much lower than the radial trap frequency, i/^ « 31 Hz. 
This behavior is quite different from that of single-species 
BEGs where excitation frequencies are always greater 
than the trap frequency [isl . [23j . and further empha- 
sizes the essentially two-component nature of the exci- 
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FIG. 2: (color) Top panels: 3D renderings of the density dis- 
tributions in a binary BEG computed from the model at 30 
and 80 ms. Each component is depicted by a contour slice 
at about half of its corresponding maximal density. Red and 
green surfaces correspond to components |1) and |2) respec- 
tively. The bottom (side) projection corresponds to the z- 
(x-) integrated density for the |1) component in our model. 
Bottom panels: radial profiles of the density distributions at 
30 and 80 ms. Symbols (squares for |1) and circles for |2)) 
correspond to the experimental data, which are averaged over 
the azimuthal angle. Solid lines depict the model results. The 
middle curve depicts half of the total density, which is nearly 
preserved during the time evolution Q]. 



tation [l8|. We also observe that the two species oscil- 
late out of phase so that the total density remains nearly 
constant in time Q . A second feature also emerges from 
inspection of the ring patterns: at 30 ms. the distribution 
of |1) atoms has a central peak ("bull's eye") which is less 
pronounced in the recurrence at 120 ms. As the system 
evolves it does not precisely revisit its earlier configura- 
tion due to atomic losses and incommensurate radial and 
axial trap frequencies. 

We used radial moments (such as 1\ = J r\ip\'^dV) to 
examine the consistency of the best theoretical values 
of the scattering lengths [131 by comparing numerically 
computed and experimentally measured density profiles. 
Through a detailed analysis of the periodicity of Ii{t), 
which represents the collective excitation leading to the 
recurrence of the ring waveform in Fig. [1] we are able to 
show that the value a22 = 95.00ao agrees well with the 
experimental data. By using the previously computed 
values ail = 100.40ao and ai2 = 97.66ao [27|, we have 
obtained remarkable agreement between the experimen- 
tal and the numerical results in all the cases considered 
up to 130 ms. At longer times, effects of the imperfect 
cylindrical symmetry of the trap become increasingly ap- 
parent in the relative position of the two condensates and 
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FIG. 3: Top panels: Experimental and numerical density pro- 
files of a 50:50 mix of |1) and |2) atoms at 30 ms for different 
A*' as indicated (in thousands of atoms). Different experi- 
mental scalings [2l[ are used to avoid saturating the images. 
Bottom panels: same as the preceding for different relative 
numbers of atoms as indicated by the |1):|2) ratio. The atom 
number is fixed at iV = 3.50(5) x lO'^, as in Fig.[Tl 



the crispness of the recurring ring patterns. 

Our model can be used to infer the complex spatial 
configurations arising from the nonlinear interactions be- 
tween the two components. We show in Fig. [2] a 3D ren- 
dering of the spatial distribution of the densities for both 
species based upon the numerical results, which reveals 
interesting aspects of the interaction effects such as the 
location of the rings in the z-direction [2^. The mod- 
eled density distributions provide a richer view of the 
separation phenomena and are more subtle than might 
be suspected from an analysis of the experimental data, 
which consist only of integrated planar projections of the 
atomic density distributions. 

We have also explored the dependence of the period 
and shape of the ring-like structures on the total initial 
number of atoms and the relative atom fraction in each 
of the two spin states. Our observations at t = 30 ms as 
we vary TV and the relative number are shown in Fig. [3l 
We find that although the precise density distribution at 
any particular time depends sensitively on both of these 
variables, the periodicity of the patterns and the quality 
of agreement between the experimentally obtained and 
the numerically computed profiles do not. 

The observed ring patterns are reminiscent of radial 
excitations, and the further we drive the system out 
of equilibrium the more rings are formed. In the top 
panels of Fig. [31 the number of rings increases as the 
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mean-field (nonlinear) interaction energy of the system 
increases. Similarly, in the bottom panels of Fig. [3l the 
number of rings increases as more atoms are transferred 
from |1) to 1 2). We understand this behavior by not- 
ing that the excess energy after transfer, as determined 
by differences in scattering lengths and trapping poten- 
tials, places the system far from its ground-state equilib- 
rium and consequently induces a long oscillatory journey 
through the system's configuration space. Stronger in- 
teractions, larger transfer ratios, and greater differences 
in the trapping potentials should therefore produce more 
dramatic dynamics. The experiment demonstrates that 
the first two conditions lead to the most pronounced ring- 
like structures; the last will be explored elsewhere. 

In conclusion, we have observed for the first time a 
largely undamped two-component collective excitation in 
which each component forms striking ring-like patterns 
that recur on a timcscalc much longer than the radial 
period. Numerical solutions to a theoretical model based 
on a modified pair of coupled Gross-Pitaevskii equations, 
along with careful experimental characterization of loss 
rates and trapping potentials, accurately and quantita- 
tively describe the dynamics of the system, and help to 
visualize aspects of the complicated component separa- 
tion behavior that are difficult to observe directly. We 
anticipate the possibility of assuming the accuracy of the 
model and using experimental data to extract scatter- 
ing length ratios from the periodic behavior of the den- 
sity distributions, much as in Ref. (20| . This technique 
would be sensitive to the interspecies scattering length, 
of crucial importance to the question of the miscibility 
of two atomic species that depends on the value of the 
determinant D = 011022 — 012 [13, U, ^ fur- 

ther extension involves self-consistently including beyond 
mean-field effects, such as thermal excitations, in order 
to match the corresponding component separation exper- 
iments at finite temperature. These directions are sub- 
jects of future work. 
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